function plot_model(model)
Mu = model.mu';
Sigma = model.Sigma;
[k,d] = size(Mu);
if d == 3
    for n = 1:k
        P = zeros(2,d);
        [V,D] = eigs(Sigma(:,:,n)'*Sigma(:,:,n));
        P(1,:) = Mu(n,:);
        for m = 1:d
            hold on;
            P(2,:) = P(1,:) + sqrt(D(m,m))*V(:,m)';
            plot3(P(:,1)',P(:,2)',P(:,3)');
        end
    end
end